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Abstract. Recent analyses of the intergalactic UV background by means of the Hen Lya forest assume that Hen and Hi 
absorption features have the same line widths. We omit this assumption to investigate possible effects of thermal line broadening 
on the inferred He ii/H i ratio 77 and to explore the potential of intergalactic He n observations to constrain the thermal state of 
the intergalactic medium. Deriving a simple relation between the column density and the temperature of an absorber based 
on the temperature-density relation T = To{l + Sy^^ we develop a procedure to fit Tq, 7, and 7/ simultaneously by modeling 
the observed spectra with Doppler profiles. In an alternative approach the temperature T of an absorber, the Hen/Hi ratio 77, 
and the redshift scale of 77 variations are estimated simultaneously by optimizing the Doppler parameters of He n. Testing our 
procedure with artificial data shows that well-constrained results can be obtained only if the signal-to-noise ratio in the He n 
forest is S/N > 20. Additionally, ambiguities in the line profile decomposition may result in significant systematic errors. 
Thus, it is impossible to give an estimate of the temperature-density relation with the Hen data available at present (S/N ~ 5). 
However, we find that only 45 % of the lines in our sample favor turbulent line widths while the remaining lines are probably 
affected by thermal broadening. Furthermore, the inferred 77 values are on average about 0.05 dex larger if a thermal component 
is taken into account, and their distribution is 46 % narrower in comparison to a purely turbulent fit. Therefore, variations of 
77 on a 10 % level may be related to the presence of thermal line broadening. The apparent correlation between the strength 
of the H I absorption and the 77 value, which has been found in former studies assuming turbulent line broadening, essentially 
disappears if thermal broadening is taken into account. In the redshift range 2.58 < z < 2.74 towards the quasars HE 2347-4342 
and HS 1700+6416 we obtain 77 a; 100 and slightly larger In the same redshift range the far-UV spectrum of HS 1700+6416 is 
best and we estimate a mean value of log 77 = 2. 1 1 ± 0.32 taking into account combined thermal and turbulent broadening. 

Key words, cosmology: observations - quasars: absorption lines - quasars: individual: HE 2347-4342, HS 1700+6416 



>> 1. Introduction 



The temperature of the intergalactic medium (IGM) preserves 
a memory of the r eionization history of the universe (e.g. 



iTheuns etaL 2002ah and p rovides constraints on th e nature of 
the reionizing sources (e.g.l Tittlev & Meiksinll2006 ). The usual 
approach to measure the temperature considers th e Doppler pa 



ramet e r distr i bution of Lyg forest absorbers (e.g.lSchaye et al 



2000[ [b 
see I The 



laye et al. 

skTEooof 



Bryan & Machacek 

(2002b) for 



199?, '200d; iRicotti et_aL 

.McDonald et al. 2001, but see iTheuns et al. (ZUUZb) tor a 
different approach). Generally, there is a power law depen- 
dence between the column density A^hi and the lowest Doppler 
parameter at this column density, i.e. the cut-off b^i„(NHi)- 
Interpreting b^in as purely thermal, the temperature of the 
IGM can be estimated. In practice, hydrodynamic simula- 
tions are required to transform the fc min-A^Hi depe ndency 
into a temperature-density relation (e.g. ISchave et aLlll999l) . 
However, even for the narrowest lines the thermal component is 
not necessarily the only contribution to the line width. Thus, the 
cut-off provides just an upper limit of the temperature. Further 
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broadening is introduced by the velocity field including pecu- 
liar differential motions and the Hubble flow. The observed dis- 
tribution of ^-parameters is also related to the amplitude of the 
primordial density fluctuations as shown by iHui & Rutledge 
( 1999) using linear perturbation theory. 

The approach of estimating the cut-off of the b(NHi) distri- 
bution to measure the temperature of the IGM is based on the 
concept of a well-define d power law rel ation between temper- 
ature and density. Hui & GnedinI (Il997i) showed this relation 
to be valid in the low-density Lya forest for uniform reioniza- 
tion models as a consequence of the interplay of photoioniza- 
tion heating and adiab atic cooling due to the cosmic expansion. 
Valageas et al.l (120021) argued that the power law relation is still 
valid even if the ionizing UV background is inhomogeneous 
as long as the gas is in local ionization equilibrium. The re- 
lation breaks at high densities since in the non-linear regime 
shock-heating due to gravitational processes becomes domi- 
nant introducing a second IGM component, the warm-hot inter- 
galactic medium (WHI M) containing the bulk of the baryons at 
the present epo ch (e.g. Cen & Ostrikei" 199^1 iDave et al ] l200lf 



Valageas et d] [2002: Kang et al.. .20051) . On the oflier hand, 



2 



C. Fechner and D. Reimers: The Hell hya forest and the thermal state of the IGM 



high-density gas may cool down in less than a Hubble time 
which is a precondition to enable star formation. Thus, cool 
dense gas exists at z ~ 2 - 3 and even at lower redshifts (e.g. 
Dave et al. Il999t IValageas et al.ll2002l) . As a consequence the 
temperature-density relation breaks at high densities (S ~ 10 
or A^Hi ~ lO'^cm"^, respectively). Above this threshold the 
gas temperature decreases with den sity (for a summ ary of the 
phases of the IGM see e.g. Fig. 1 of lValageas et al ]^002). 

In agreement with the theoretical scenario iMisawa et al 



(|2004 observed, indeed, that due to shielding effects high den- 
sity lines are narro wer and therefore c older than low density 
absorbers. Whereas iRauch et al.l (l2005h found from an inves- 
tigation of the velocity structure of Lya forest clouds in the 
spectra of multiple QSO sight lines that the large scale mo- 
ti ons in the IGM are d ominated by the Hubble flow. According 
to iRauch et al.l (120051) some absorbers undergo early stages of 



gravitational collapse but the majority of absorbers apparently 
arise from expanding filaments accelerated by the infall into 
more massive structures. 



In a study of the reionization of helium lGleser et al.l (120051) 
find that the temperature-density relation may trace the stages 
of He n reionization. With ongoing reionization the slope of the 
temperature-density relation gets steeper and a significant scat- 
ter is introduced. At z ~ 2-3 the scatter affects in particular the 
low density regime. When helium reionization is completed the 
slope of the power law flattens again and the amount of scatter 
decreases since adiabatic cooling dominates over photoioniza 



tion he ating. Similar results have been obtained bv lBolton et al 



(12004 who examined the effects of radiative transfer on the 
temperature-density relation . Indeed, the mean temperature of 
the I GM peaks at z ~ 3 dSchave et al.1 12OOOI: iTheuns et af 



l2002aiibi) which is thought to be due to additional heating be- 
cause of He II reionization. 

From the analysis of the resolved He 11 Lya forest in the 
spectra of the quasars HE 2347-4342 and HS 1700+6416 a cor- 
relation between the density of the absorbers and the hardness 
of the ionizing radiation (represented by the parameter 77 which 
is the He ii/H i column density ratio) has been found in the sense 
that higher column densit y absorbers are apparently exposed to 



a harder radiation field (iShull et al.ll2004t iZheng et al ■112004 



Reimers et al. 2006; Fechner et al. 2006a). There are several 



suggestions for a physical explanatio n of this correlation like 



self-shielding and shado wing effects (IMaselh & Ferrarall2005 



Tittlev & Meiksinll200^ and the spatial distribution of QSOs 
as the sources of ha r d radi ation ( Bolton et al. 2006) . However, 
Fechner & ReimersI ( 20061) suggested that part of the correla- 
tion might be due to the unjustified assumption that the line 
widths of the absorption features are dominated by purely tur- 
bulent broadening. If this assumption is incorrect, a correlation 
between the 77 value and the absorber density will be introduced 
artificially. In order to reproduce a given equivalent width of an 
absorption line the column density will be overestimated if the 
line width is underestimated. This effect is in particular pro- 
nounced for strong absorption features with high column den- 
sities that are beyond the linear part of the curve of growth and 
thus introduce an apparent correlation between 77 and H i col- 
umn density. 



Measurements of observed f eatures indicate th at the lines 
are indeed turbulently broadened dZheng et al .|2004 . However, 
this result is based on a limited statistical sample of 8 absorbers. 
Liu et al 1 ([2006) claimed to confirm this finding using hydro- 
dynam ic simulations. On the other hand, iFechner & Reimers 
(120061) showed that statistically, lines with apparent low rj val- 
ues are often dominated by thermal broadening. In order to 
solve this apparent conflict further investigation is needed. 

In this work we introduce a simple approach to consider the 
temperature of the IGM when analyzing the He n Lya forest. 
This provides, at least in principle, an independent method to 
measure the temperature-density relation (often referred to as 
equation of state of the intergalactic medium). Since HE 2347- 
4342 and HS 1700-1-6416 are the only objects where high reso- 



lution FUSE spectra are available ([K riss et al. 2001t lShull et al 



2004t IZheng et al.ll2004i iFechner et al.. .2006a,). we apply our 
method to these sight lines. However, as we will see, the data 
quality (S /N ~ 5) is insufficient to derive tight constraints for 
the temperature-density relation. 

After shortly reviewing several equations important in 
the context of the temperature-density relation in Section |2] 
Section [3] presents the methodical approach to derive the pa- 
rameters of the temperature-density relation and 77 from ob- 
servations of the Hen Lya forest in comparison to Hi. This 
method is tested using artificial data in Section |4] and applied 
to the observations in Section |5] In Section |6] we introduce an 
alternative approach to estimate the temperature for the indi- 
vidual absorption lines more directly and discuss the effect on 
the inferred 77 values. 

2. The thermal state of the IGM 

Hui & GnedinI (1997) have shown that there is a well-defined 



power law relation between temperature and density, often re- 
ferred to as equation of state for the low density IGM {S < 10, 
i.e. A^Hi < lO'^cm-^). 



T = Tod +6) 



(1) 



where 5 = | - 1 is the gas overdensity and Tq the temperature 
at the mean density. This equation of state can be converted 
into a power law relation between the cut-off Doppler param- 
eter fcmi n and the H i colum n density using cosmological simu- 
lations ( Schaye et al.lll999 ). Therefore, estimates of the param- 
eters To and 7 are usually based on measuring the cut-off of 
the b{NHi) distribution. The combination of Hi and He 11 data 
would provide, at least in principle, an independent estimate 
of the thermal state of the IGM. In order to relate the equation 
of state to the observed H i and He n spectra it is necessary to 
connect the H i column density A^h i and the temperature of the 
absorbing material. 

Assuming that the almost fully ionized gas is in photoion- 
ization equilibrium, the number density of neutral hydrogen 
can be written as 



72HI * 1.41 -lO-'^cm-^d +5)2(1 +z)^ (y^) 



-0.755 



-I 



0.02 



9x2 



(2) 
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with the photoionization rate of hydrogen Fhi- In order to de- 
rive this relation we have used 



nil = 



(i-y)(i+z)^(i + 5), 



(3) 



with the primordi al mass fraction of helium Y ^ 0.244 
jBurles et al.ll200lb . 

Since A^hi = ■ I, Eq. Q leads to a relation between Hi 
column density and overde nsity 6. However, the length / of the 
absorbers has to be known. iBrvan & Machacekl (l2000h choose 
a typical size of absorbe rs from simulations scaling with I oc 
(1 +z)-'.(Sd^ (E()Oll) argued that the characteristic size of 
an absorber along any line of sight will generally be in order of 
the local Jeans length 



(4) 



with the fraction fg of mass in the gas which should be close 
to the universal value fg ^ Qb/Qm ~ 0.16 in the low density 
Lya forest. Combining Eqs. Q, and (|4| leads to an analytic 
expression for the relation between the observable column den- 
sity of ne utral hydroge n A^h i and the corresponding overdensity 
(see also ISchavell200lh 



A^Hi ^ 6.50-10"'cm--(l+5)'-5(l+z)^-^(^) 



-0.255 



Thi 



-1 /r> ;2\l-5 



0.02 



0.16 



0.5 



(5) 



The latter three parameters are in the order of unity. Then, in- 
serting the equation of state yields 

logA^Hi ~ 11.833 + 4.5 ■log(l +z)- log To 

r- 1 



1.5 
7^ 



0.255 ■ log r , 



(6) 



where A'hi is measured in cm"^ and the temperatures are given 
in K. Thus, we have found a relation between the H i column 
density and the temperature of an absorber which only depends 
on the redshift and the parameters of the equation of state y and 
log Tq. However, its validity is limited to gas with low and mod- 
erate densities. At higher densities (A^hi S lO'^ cm"^) radiative 
cooUng dominates over adiabatic cooling due to the Hubble 
expansion, leading to lower temperatures. Furthermore, a pos- 
sible scatter in the temperature -density relation due to helium 
reionization ( Gleser et al JI2005 ) is neglected. 



3. Method 

Three different analysis methods have been applied to the ob 
servations of the Hen Lya forest so far: line profile fit ting 
dKriss et al. l l200ll: IZheng et alJl2004l: iFechner et alj|2006ah. an 



apparent optical depth method (Shull et al.'2004' 'Fechneretal 
[2006a) and spectrum fitting (Fechner & Reime rs 2006). The 
latter fits the high-quality optical H i spectrum directly to the 
He II data, and allows to estimate the length scales of rj varia- 
tions. Among these three methods only profile fitting provides 



the possibility of taking into account different line widths of H i 
and He ii as would be introduced by thermal line broadening. In 
order to preserve the advantage of the spectrum fitting method 
to estimate scales on which the ionizing background fluctuates, 
we combine the traditional profile fit analysis and the spectrum 
fitting method. 

The H I Lya forest is fitted with Doppler profiles. From the 
resulting line list an artificial He ii spectrum is computed and 
convolved with an appropriate instrumental profile. The arti- 
ficial spectrum can then be fitted to the FUSE data by scal- 
ing the H I column density. Furthermore, the line width of the 
He II features can be addressed via the Doppler parameter /^Hen 
which could lead in principle to an independent estimate of 
the temperature of the IGM and the parameters of the equation 
of state. As we will see below, fitting the temperature is sen- 
sitive to the noise level of the Hen spectrum. The quality of 
present day He ii data from FUSE is insufficient to derive well 
constrained temperatures and scales of variation for rj simulta- 
neously. Consequently, we will a priori select redshift ranges 
where rj is assumed to be constant. 

Fitting line profiles introduces fu rther uncertainties. As 
demonstrated bv lFechner et al.l (l2006al) . a scatter of 10 - 15 % 
in T] is expected even if the underlying value is constant. This is 
due to the fact that any line profile model recovers the observed 
data only imperfectly. Additionally, the high noise le vel of the 



FUSE data introduces further uncertainties (see also iLiu et al 



|2006). Since these problems are inevitable, we expect that our 
results can be interpreted statistically but conclusions for single 
absorbers should be considered very carefully. 

The modified spectrum fitting method works as follows: 
First, the H i data is fitted and a line list with the values of red- 
shift z, column density log A^hi, and Doppler parameter bni is 
prepared. Then values of the parameters log Tq and y are se- 
lected. The temperature T of each single line is computed ac- 
cording to Eq. (|6]l. The Doppler parameter of He ii is then given 
by 



^Heii = Jbl^^-2kT-{- ^1, 



(7) 



where k is Boltzmann's constant and mu and mee are the 
masses of the hydrogen or helium atom, respectively. If the 
temperature is very high, the term beneath the square-root in 
Eq. d?) may become negative. In this case the Doppler param- 
eter of H I is assumed to be completely thermal and the tem- 
perature is set to r - /j^j ■ niYililk). From the adjusted line 
parameters an artificial Hen spectrum is generated and fitted 
to the dat a using a -procedure where log rj is varied (as de- 
scribed in .Fechner & Reimers 2006) . This means, we assume 77 
to be constant in the considered redshift range. This procedure 
is repeated for several values of log Tq and y. We explore the 
parameter space 3.50 < log To < 4.50 and 0.70 < y < 1.80 
where realistic results are expected. By comparing the result- 
ing ;t'^-values for each combination of log Tq and y the best fit 
parameters can be found. In total three free parameters (log rf, 
log Tq, and y) are estimated for a given redshift interval. 

Eq. O assumes that the thermal and turbulent Doppler pa- 
rameters both lead to a Gaussian profile and combine to the 
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total Doppler parameter as - b^^^^^ + b^^^^- Since the tur- 
bulent part of the ^-parameter summarizes all kinds of veloc- 
ities except thermal, the broadening profile is not necessarily 
Gaussian. In particular, due to the atomic parameters of helium, 
absorbers with moderate He ii column densities probe the low 
density Lya forest and are probably dominated by the differ- 
ential Hubble flow rather than microscopic turbulence. Thus, 
Eq. (|7]i represents the most simple approach to combine ther- 
mal and turbulent line broadening. This should be sufficient for 
a first analysis. 

4. Application to artificial data 

In order to explore the potential of the introduced approach 
we examine artificial data which is created on the basis of 
the statistical properties of the Lya forest. The column den- 
sity distribution functio n with /3 - 1.5 is adopted from 
Kirkman & Tvtleild 19971) . Values in the range 1 1 .0 < log A^h i < 



18.0 are simulated. The Doppler parameter distribution is de- 
scribed by a truncated Gaussian (bmin - lOkms"') centered 
at 27kms~' w ith a w ith of cri, - 8.75 km s"' in agreement 



with Huetal 



(Il995h . In addition the temperature is linked 
to the Hi column density according to Eq. The param- 
eters of the temperature-density relation are se t to the real- 
istic values log To = 4.3 and y = 1.3 (e.g. ISchaye et al 
Furthermore, w e chose log// - 1.903 (i.e. rj 



e-g- 



80; 



Kriss et al 



2001h . The spectra cover the redshift range 
2.56 < z < 2.92 with R = 40000 in the Hi and = 15 000 
in the He ii Lya forest. In order to study the importance of the 
data quality several signal-to-noise ratios are chosen for the ar- 
tificial Hen spectrum (S /N = 5, 10, 20, 30, 50, 100). The 
for the Hi Lya forest is set to 100 which is comparable to the 
noise level of the best optical spectra available today and is 
sufficient to derive a reasonable line list from the fit. 

The H I Lya forest is fitted with Doppler profiles once and 
afterwards all artificial He ii spectra are fitted on the basis of the 
derived line list according to the method presented in Section 
[3] Fig. [T] shows the 1 cr and 2 cr confidence contours for the pa- 
rameters log T() and y of the temperature-density relation (Eq. 
[T]). The point of the best fit is marked by a star, while the true 
underlying value is indicated by a cross. As expected, the re- 
sults are better constrained if the S /N of the He ii spectrum is 
increased. In case of S /N - 5 a reasonable estimate cannot 
be derived since the variation of the absorber's temperature af- 
fects the wings of the line profiles only marginally and due to 
the high noise level very similar -values are produced. In or- 
der to derive reasonable constraints for the temperature-density 
relation a signal-to-noise ratio of at least 5 /A^ ~ 20 ... 30 is re- 
quired for the Hen data. The results are summarized in Fig. 
|2] which presents the best fit values and their 1 cr uncertainties 
in comparison to the true underlying values. It becomes clear 
that also for S /N - 100 the derived log Tq and y are underes- 
timated. The reason for this systematic effect is the imperfect- 
ness of the adopted line list (see below). However, the inferred 
J] values are in good agreement with the underlying rj = 80. 

For comparison we also perform a fit assuming pure tur- 
bulent line widths, i.e. buen - ^hi- The resulting rj values are 
presented as triangles in the right panel of Fig. |2] Obviously, 



4.2 
4,0 
3. a 



S/N = 5- 
" " 


S/N = 10- 
~ 


S/N = 20- 
■ N " 


1 1 1 i_ 

"s/N = 30- 
" "'i. " 


_i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i_ 
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" " 


_i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i_ 
S/N = 100- 



Fig. 1. Contour plots for the parameters of the temperature- 
density relation derived from the simulated data. The crosses 
represent the position of the true values of y = 1.3 and log To - 
4.3 while the stars indicate the estimated minimum. The outer 
contours give the joint 1 cr (solid line) and 2 cr (dotted line) con- 
fidence levels, respectively. The inner contour represents the 
1 cr confidence level for one free parameter, i.e. the projection 
of the inner contour to the axes leads to the 1 cr errors of log Tq 
and y, respectively. The signal-to-noise ratio of the artificial 
He II spectra is indicated in each panel. The confidence levels 
in case of S /N - 5 (upper left panel) are outside the presented 
parameter range. 



i: 



30 40 60 BO 100 



20 40 60 ao 100 



Fig. 2. Results for the fit of the artificial data. The panels show 
the derived values for log To, y, and log rj (from the left to the 
right). The error bars correspond to the inner 1 cr confidence 
contours of Fig. [T] The dotted lines indicate the true values. 
For comparison the 77 values estimated under the assumption of 
pure turbulent line broadening are presented as triangles in the 
right panel. 



the rj value is underestimated if pure turbulent line broaden- 
ing is assumed. This finding is essentially independent of the 
signal-to-noise ratio of the He 11 and also of the modifications 
of the line list (see below). On average the pure turbulent 77 
value is ~ 0.05 dex lower than the 77 derived when considering a 
temperature-density relation and thus about 0.05 dex lower than 
the underlying value. This is still consistent within the 1 cr level 
in case of low S/N but the underestimation will be significant 
as the error bars decrease for higher quality data (S /N > 50). 

As noted above, the best fit depends on the quality of the 
adopted line list. This is illustrated in Fig. [3] We performed 
fits to the artificial Hen spectrum with S/N = 100 based on 
three different Hi line lists. The "perfect" line list includes 
all absorption fines that have been generated to compute the 
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Fig. 3. Contour plots for the parameters of the temperature- 
density relation derived from the artificial data using three dif- 
ferent Hi line lists for the Hen fitting (see text). The crosses 
represent the position of the true values of y = 1.3 and log Tq = 
4.3 while the stars indicate the estimated minimum. Contour 
levels are given as in Fig. [T] 



artificial spectra. In this case the correct parameters are, of 
course, reproducible (Fig. [3^). We find log Tq - 4.30 + 0.02 
and J - 1 .30^y yg. For a "reasonable" line Ust based on a thor- 
ough fit of the H I spectrum we obtain log Tq - 4.25 + 0.05 and 
y - 1 .24 + 0.04 (Fig.[3j'). In particular y is not reproduced cor- 
rectly although the estimated value has apparently small error 
bars. In this case the imperfectness of the line list introduces 
additional systematic uncertainties. 

Furthermore, we examine a realistic but "problematic" case 
(Fig. [3};). The difference between the realistic and the problem- 
atic line list is illustrated in Fig.|4] A broad, saturated H i feature 
has been fitted with two Doppler profiles (realistic case) or with 
one component only (problematic case). The upper panel of 
Fig.|4]shows only very slight discrepancies between the result- 
ing profiles. However, there are more prominent differences in 
the He ii fit (lower panel of Fig. |4]i. As a consequence of the ab- 
sence of a second component the feature has to be modeled by 
a broader profile to achieve a reasonable fit. Therefore, the tem- 
perature T has to be lower according to Eq. (|7]i. Because of the 
given H i column density of the line y has to decrease (see Eq. 
ISll. Indeed, we find log To = 4.25!° °^ and y = 1.08!° °' (Fig. 
[3};). Additionally, a weak H I line (logA^Hi - 1 1 .7) leading to a 
more prominent Hen feature (logA^Hen = 13.4) at 1001.0 A is 
missed in both models introducing further deviations from the 
optimum model. 

However, in case of observed spectra the best strategy to 
produce the H i line list would also take into account higher or- 
der features of the Lyman series. Thus, the estimation of the 
line parameters of a saturated absorption complex would be 
more accurate than in this simple simulated case. Due to blend- 
ing or incomplete spectral coverage a few features with am- 
biguous line parameters might still be left unsettled even for ob- 
served spectra. But the existence of unidentified, weak H i fea- 
tures below the detection limit can never be excluded. For this 
reason profile fit based analyses of the He ii Lyo' forest include 
He n lines without detected H i counterp arts (Kriss et al. .200 1; 
Zhenp et al. 12004 Fechner et al.ll2006al) . One should keep in 



mind that even one missing component may affect the results 
significantly. 
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Fig. 4. Illustration of the importance of the H i fit. The upper 
panel shows a portion of the artificial H i data (histogram-like) 
and two different fits where the saturated feature at 4007 A is 
modelled with one (dashed line) or two (solid line) Doppler 
components, respectively. The inset in the upper panel illus- 
trates the only marginal difference between both fits. The lower 
panel shows the artificial Hen data (histogram-like) and the 
best fits in case of the one component approach (dashed line; 
log To ~ 4.25, y ~ 1.08) or two component fit (solid line; 
log To ~ 4.25, y ~ 1.24), respectively. The dot-dashed line 
represents the best fit using the perfect line list (log Tq - 4.30, 
y = 1.30). 



5. Application to observations 

Though the data quality for the two He n quasars observed with 
FUSE, HE 2347-4342 and HS 1700H-6416 (S/N ~ 5), is well 
below the required noise level of S /N > 20, we apply the 
method presented in Sect.|3]to the observed spectra. Our main 
intention is to explore whether any conclusions can be inferred 
even from this noisy data. In order to avoid additional degrees 
of freedom the scales to perform a fit are selected a priori. We 
select spectral ranges with a size of ~ 5 A, corresponding to 
Az ~ 0.016. Scales of variation of the U V background in this 



order o f magnitu de have been estima ted b ylFechner & Reimers 



(120061) . whereas IShuU et al.1 (120041) and iBolton et al 



(l2006h 



find much smaller scales in the range Az ~ 0.001. Realistically, 
we do not expect to be able to derive robust constraints on the 
parameters of the temperature-density relation. However, this 
approach should be sufficient to roughly estimate the effect of 
partly thermal Doppler parameters on the results. 

The far-UV d ata of HE 2347-4342 are presented in 
Zheng et alJ (120041) . We use the normalization described in 
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Fig. 5. Results from HE 2347-4342 for log 77 (open circles; up- 
per panel), log To (middle panel), and y (lower panel). The tri- 
angles in the upper panel present the resulting rj value assuming 
pure turbulent broadening. For clarity error bars are not given 
for the turbulent case, but are of the same size as those of the 
temperature-density relation fit. The 1 cr uncertainties of log To 
and -y would comprise the whole presented parameter range. Fit 
sections of ~ 5 A have been selected a priori. Sections marked 
with a cross represent redshift intervals without extremely sat- 
urated absorption features. The dotted lines in the lower pan- 
els indicated the expected value interpolated from the measure- 
ments of .Schaye et al.. (.2000.) . 



Fig. 6. As Fig.|5]but forHS 1700H-6416. 



Fechner & Reimerll ( 2006 ). The spectrum is binned to 0.05 A 
which corresponds to the resolution of R - 20000 at 1000 A 
(and the actual value ~ 15kms ') and reduces the noise to 
a S /N ~ 5 level per bin in the whole spectrum. In order to 
avoid further uncertainties we exclude the spectral range con- 
taining Ly/3 or higher order Lyman series lines. For the Hi 
Lyg forest towards HE 23 47-4342 we adopt the line list of 
Fechner & Reimers (l2006h based on a Doppler profile fit of 
high-quahty VLT/UVES spectra (R ^ 45 000, S /N ~ 100). We 
slightly modify the parameters of the Lyman limit system com- 
plex at z ~ 2.735 which are now estimated using the Lyman 
series up to Ly-9. 

The FUSE data of HS 1700-H6416 are presented in 
Fechner et aP (l2006 al). It has a resolution of R ~ 20000. The 
signal-to-noise ratio is ~ 5 and slightly better. The line list 
for H I Lya forest is adopted from the results of Fechner et al.l 
(12006a') who use high-quality data from Keck/HIRES (R a; 
38 500, S /N ~ 100). Furthermore, we correct for metal lines 



in the FUSE spectral range as described bylFechner & Reimers 
(12006') using the models of ' Fechner etalT (l2006b'). 



As a compromise between constructing a model as reahstic 
as possible and the caveats due to the limited quality of the data 
we fit the selected spectral ranges assuming constant values of 
77, log To, and 7. For comparison we also estimate a fit assum- 
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Fig. 7. Correlation between the difference in the quality of the 
fit, i.e. the;^'^, and the difference in the rj value. The better the fit 
taking into account a temperature-density relation in compari- 
son the the fit assuming pure turbulent broadening, the larger 
is the resulting log 77 for the best fit. The dotted line repre- 
sents a linear fit to the Lye forest data, i.e. HS 1700+6416 and 
HE 2347-4342 at redshifts z < 2.75. 



ing pure turbulent line widths. The results are presented in Figs. 
|5](HE 2347-4342) and|6](HS 1700-H6416). Towards HE 2347- 
4342 the fit inte rvals at z > 2.75 c over the range of patchy 
He II absorption (iReimers et al ] ll997l) which marks the tail end 



of the epoch of H e 11 reionization (also e.g. iKriss et alJl2001 



Zheng et alj|200 4'). For these redshift ranges fur ther uncertain 



ties like the incomplete reionization of He 11 (e.g. lReimers et al 



19971 l2005h or simply the fact that the FUSE spectrum shows 
completely saturated, broad absorption troughs, lead to ex- 
tremely unstable fits which should not be considered further. 
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Fig. 9. Distribution of the parameters of the temperature- 
density relation estimated individually for the sight lines to- 
wards HS 1700+6416 (black) and HE 2347-4342 (light; z < 
2.75). The value of log To - 4.3 has been fixed when estimating 
y (left panel), while y = 1 .3 when log To has been optimized 
(right panel). 



Generally, the r] values in the pure turbulent cases are 
slightly lower than if the thermal state of the IGM is taken into 
account, consistent with the results from the artificial data. On 
average the turbulent 77 is lower by roughly 0.05 dex consider- 
ing the Lya forest only (i.e. z < 2.15). The difference in the rj 
value is correlated with the quality of the fit, i.e. with the 
which is shown in Fig. |7] If a thermal component to the line 
width is considered, the He 11 lines become narrower in com- 
parison to the pure turbulent case. In order to produce roughly 
the same equivalent width the column density of He 11 and thus 
77 has to increase. The correlation presented in Fig. Qmeans if 
considering the thermal state of the absorbers significantly im- 
proves the fit in a selected He 11 wavelength range, the inferred 77 
value will be larger. And even if the fit is improved only slightly 
^mrb ^ ^ value increases on average by 0.02 dex 

(median 0.01 dex). However, the deviations are still within the 
1 cr uncertainties of the estimate (see also the upper panels of 
Figs. |5] and |6]l. 

The parameters of the temperature-density relation esti- 
mated by the fits are rather unconstrained. As described in Sect. 
[3] we explored the parameter space 3.50 < log < 4.50 and 
0.70 < y < 1.80. Within these limits all values produce models 
fitting to the data on a 1 cr level. Therefore, no error bars are 
presented in Figs. |5] and |6] Nevertheless, some tendencies are 
noteworthy. 

The vast majority of the fit intervals exhibit an extreme 
value for y (< 0.7 or > 1.8) and/or log To (< 3.5 or > 4.5). 
Excluding the redshift range with patchy He 11 absorption (z > 
2.75) from further discussion, 81 % of the fitting sections yield 
extreme values of the temperature-density relation. Among 
these cases 34 % lead to a temperature-density relation with 
y> 1.8 and log To < 3.5. 

In order to give an interpretation of the temperature-density 
relation in this case some examples are illustrated in Fig. [8] A 
large value of y leads to a steep relation between column den- 
sity and temperature. Consequently, the temperature of high 
column density absorbers is much higher than the tempera- 
ture of low column density absorbers. The parameter log To 



rules the absolute values of the temperature. The lower log Tq, 
the lower is the general temperature level. Finding y > 1.8 
means that the line widths bneu of high column density ab- 
sorbers are significantly narrower than bui- Therefore thermal 
line bro adening is favored. This ide a may support the conclu- 
sions of iFechner & Reimers (l2006l) who suggest that thermal 
line broadening might be important for high column density 
absorbers. 

For about 9 % of the extreme cases we find y < 0.70 and 
log To < 3.50, the lowest values in the parameter range ex- 
plored. At this point the lowest temperatures are produced for 
high column density absorbers since the slope (y - 1) of the 
temperature-density relation gets negative and log To -3.5 sets 
the absolute scale to low temperatures. Thus, thermal broad- 
ening might play a minor role for the line widths of the H i 
and Hen features. For these fit sections the 77 values estimated 
from the pure turbulent and the thermal state model are nearly 
identical supporting the idea that turbulent line broadening 
is dominating in these redshift ranges. At z ~ 2.43 towards 
HS 1700-1-6416 both derived 77 values are heavily overproduced 
due to problems with two saturated, broad absorption features 
that dominate the selected spectral portion. 

The results are expected to be more stable if regions with 
saturated or extremely weak H i are excluded. The data points 
marked with a cross in the lower panels of Figs. |5] and |6] 
represent the sections free from strongly saturated features. 
However, the results presented above for the whole spectra is 
also valid when considering only the selected portion of the 
data. Differences should be expected for data of higher quality. 

Furthermore, we test to constrain the parameters of the 
temperature-density relation by fixing one parameter at a rea- 
sonable value and varying the other one and log 77 simulta- 
neously. The estimated distributions of log 77 are very simi- 
lar to those presented in Figs. |5] and |6] The rms-deviation 
from the distributions shown is ~ 0.15 or even smaller Fixing 
log To = 4.3 and varying y still leads to inconclusive estimates 
of y for both lines of sight. The derived distributions are pre- 
sented in the left panel of Fig.|9] The majority of redshift bins 
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prefers an extreme value y > 1.8 like in the standard pro- 
cedure. Furthermore, the distribution might peak at y ~ 1.2, 
which would be consistent with independent estimates from 
the hterature (see below). However, the detection is insignifi- 
cant for the present data. Similar results are found if y = 1.3 
is fixed. The estimated distributions of log To are shown in the 
right panel of Fig. |9l In this case low values of log Tq are pre- 
ferred like in Figs. |5] and |6l There might be a slight tendency 
for log To ~ 4.4 considering the HE 2347-4342 line of sight, 
while towards HS 1700H-6416 log To ~ 4.1. However, these re- 
sults are rather inconclusive. In order to give more significant 
constraints. He ii data of better quality is needed or at least a 
more reliable statisti cs based on a larger sample . 

In t he literature iBrvan & Machacekl fcoOOl). iRicotti et al 



( I2OOOI) . ISchave et all i200(f). and lMcDonald et al.1 j200lh mea- 
sured the parameters of the equation of state in the consid- 
ered redshift range (2.3 < z < 2.9) estimating the lower 
cut-ofF of the b(Niii) distribution. Their resulting values are 
in the range 0.8 < y < 1.7 and 4.05 < log To < 4.45 
with a strong redshift dependence. Because of the just com- 
pleted epoch of He 11 reionization y increases steeply towards 
lower redshifts whil e log Tp decreases wi th z- Fitting a straight 
Une to the results of Schave et al. ( 2000l) in the redshift range 
1.8 < z < 3.2 we find y = -(0.37 + 0.16) ■ z + (2.15 + 0.42) 
and log T o ^ (0.27 ± 0.06) ■ z + (3.50 + 0.15), respectively. The 



results of Schaye et al. ( 2000l) are selected since they cover a 
large redshift range including the portion we are interested in. 
The other authors use slightly different methods and provide 
fewer data points (difference between the methods and system- 
atic offsets in the results are dis cussed in Brvan & Machacek 

EOOQ: Ricotti et al 2000: McDonaid"et al.ll200lh . For compari- 

son the fit to the Ischaye et al.l (l2000h data is indicated in Figs. 
|5]and|6]as dotted lines. It is obvious that our results scatter over 
the range given by the literature and beyond. Furthermore, there 
is no significant co ncentration of measurements towards the 
Schave et al. interpolation. Therefore, a detailed com- 



parison is meaningless. As demonstrated in Section|4]UV data 
of higher quaUty (S /N > 20) would be needed to derive tighter 
constraints for the temperature-density relation using the He 11 
Lya forest. 

In ad dition, recent indi c ations for a small value of y w 1.1 
found by D'Odorico et al. ( 2006 ^ investigating the clustering 
properties in a sample of close QSO pairs, or even hints of an 
inverses temperature-density relation (y * 0.5) which would 
improve a fit of the observed Lya flux probability distribution 
functio n of a sample of 63 quasars reported by Becker et al. 
(l2006h cannot be addressed here due to the limited data quality 
of the He 11 spectra. 

6. Direct temperature estimate 

Even though we showed in the previous Sections that an es- 
timate of the IGM temperature requires better UV data, we 
also follow an alternative approach to estimate the tempera- 
tures of the absorbers more directly. An advantage of this alter- 
native procedure is that the scale of the UV background vari- 
ations can be estimated. Therefore, the method described in 
Section |3] is adjusted. Instead of the temperature-density rela- 
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Fig. 10. Best fit towards HE 2347-4342 in the redshift range 
2.58 < z < 2.73 (upper panel). Fitted parameters are the tem- 
peratures of the absorbers (lower panel), the He n/H i column 
density ratio t] (middle panel), and the scale of 77 variation. The 
resulting temperature are presented as open circles in the lower 
panel. However, most of them have large error bars and are ac- 
tually upper limits. No temperature is given if turbulent line 
broadening is preferred (see also text). At z ~ 2.69 two sat- 
urated H I lines appear to be stronger than the corresponding 
He II features in the line center leading to a poor fit in this range. 



tion the temperature of each line itself is fitted by varying the 
Doppler parameter of He 11 according to Eq. (|7]i. We start with 
the pure turbulent case bueu - bui and then change the ther- 
m al c omponent. Furthermore, a scale estimation as described 
in|Fechner & Reimers (2006) is implemented. 

A fit is computed for a given size of an interval. Then the 
bin size is increased until the;^'^ of the best fit gets worse again. 
While the spectrum fit method uses the observed optical data 
directly and may therefore facilitate a pixelwise increment of 
the fit interval, the approach presented here is based on a super- 
position of Doppler profiles. Therefore it has to be ensured that 
all components possibly contributing at a certain wavelength 
are taken into account. Experimenting with several strategies 
in test calculations we choose a starting interval of 0.5 A in 
Hen which is increased by 3 pixels, i.e. 0.15 A, per iteration. 
Furthermore, if the enlargement of the fit interval leads to a 
larger the interval is increased (up to 5 times) to test whether 
a better fit can be obtained. 

The best fit for the redshift range 2.58 < z < 2.74 towards 
HE 2347-4342 is presented in Fig. [TO] The corresponding fit 
of the spectrum of HS 1700H-6416 is shown in Fig. [U] The 
redshift range 2.58 < z < 2.74 is chosen since the quality of 
the FUSE data in the corresponding spectral range is best for 
both lines of sight. In the following we will discuss the results 
for this regime. 

As can be seen from the lower panels of Figs. [TO] and 
fTTl the estimated temperature is only poorly constrained. 
Temperatures can be estimated for about 55 % of the lines (pre- 
sented as circles). However, in most of the cases (65.5 %) the 
1 cr uncertainties cover the whole range from the high temper- 
ature cut-ofF, when bui is interpreted as completely thermal, to 
the pure turbulent limit. On the other hand, 45 % of the lines 



C. Fechner and D. Reimers: The Hell hya forest and the thermal state of the IGM 




3.58 2.59 2.60 2.61 2.62 2.63 2.64 2.65 2.66 2.67 2.68 2.69 2.70 2.71 2.72 2.73 2.74 



2.58 2.59 2.60 2.61 2.62 2.63 2.64 2.65 2.66 2.67 2.6B 2.69 2.70 2.71 2.72 2.73 2.74 



III 



mm 



I If I 



2,58 2.59 2.60 2.61 2.62 2.63 2,64 2,65 2,66 2,67 2,68 2,69 2,70 2,71 2,72 2,73 2,74 

Fig. 11. Best fit and resulting parameters towards 

HS 1700+6416 in the redshift range 2.58 < z < 2.75; 
lines and symbols as in Fig.fTOl 
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Fig. 12. Distribution of the estimated parameters in the red- 
shift range 2.58 < z < 2.74 towards HE 2347-4342 and 
HS 1700-H6416. Presented are the distributions of the line tem- 
perature log T (left panel), the He n/H i column density ratio 
log 77 per scale (middle panel, thick line), and the scale size 
Az of rj variations (right panel, thick line). The thin-lined his- 
tograms in the right panels indicate the distributions of log rj 
and Az, respectively, in the case of turbulent line widths bueu = 

favor turbulent line broadening or yield a very low tempera- 
ture (we set the minimum temperature to lOK corresponding 
to a thermal line width < 0.5 km s ' which can be interpreted 
as purely turbulent). The upper limits estimated for those ab- 
sorbers are also indicated in Figs. [TOlandfTTl In total only 6 
lines (~ 2 %) lead to well-constrained temperature estimates. 
All of them favor the purely thermal interpretation of the H i 
Doppler parameter. The resulting distribution of the line tem- 
peratures is shown in the left panel of Fig. [12] excluding com- 
ponents that prefer turbulent line widths. The average value is 
log r = 4.41 +0.46. 

In order to obtain an estimate of the equation of state the 
derived temperature is plotted versus the H i column density in 
Fig. [13] Due to the large confidence intervals of the temper- 
ature only a slight correlation between the log A^h i and log T 
can be seen. The Spearman rank-order correlation coefficient 
is rs = 0.22. Consequently, a linear fit to the data accord- 
ing to Eq. (O is rather meaningless. The dotted line in Fig. 
[T3]indicates a temperature-density relation with the parameters 
y = 1.17 and log To = 4.23 as extrapolated from the results of 
ISchaye etal (i2000. see previous section). Our temperature es- 
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Fig. 13. Line temperature log T versus H i column density 
logA^Hi for the combined sample of HE 2347-4342 and 
HS 1700+6416 in the redshift range 2.58 < z < 2.74. Only 
lines that lead to realistic best-fit temperatures are indicated. 
The mean redshift of the considered line sample is z = 2.663 + 
0.082. The dotted line indicated a temperature-density relation 
w ith y = 1.17 and log T o = 4.23 as interpolated from the results 
of ISchave etaP (I2OO0I) . 



timates have huge error bars (the majority of points presented 
in Fig. [T 3] are in fact upper limits), but the relation derived 



from the ISchave et al. I (l2000l) results is not inconsistent with 
our measurements. Harder constraints of the line temperature 
are needed to improve the results. Also the direct temperature 
estimate would benefit from higher S /N data of the He n Lya 
forest. 

Another point is the determination of the column density 
ratio 77 and its variation with redshift for the pure turbulent 
broadening versus the thermal plus turbulent broadening cases. 
The right panels of Fig. [12] show the distribution of the re- 
sulting 77 values and the scales of its variation. The average 
Heii/Hi ratio for the combined sample is log 77 = 2.16 ± 0.62 
(median 2.01), where the line of sight towards HE 2347-4342 
contributes most of the scatter. Towards this QSO we find the 
mean value log 77 - 2.21 + 0.82 (median 1.88; note the outliers 
in Fig.[IO]l while HS 1700+6416 yields log 77 = 2.1 1 ±0.32 (me- 
dian 2.03). Also shown in Fig.[T2]is the distribution of 77 values 
for the case of fixed /jhcii = ^hi Doppler parameters. The pure 
turbulent fits yield roughly 0.1 - 0.2 dex lower 77 values. For 
the combined sample we find the mean log 77 = 1.98 ± 0.59 
(median 1.90). The results for the individual sight lines are 
log 77 = 2.00 ± 0.63 (median 1.85) towards HE 2347-4342 and 
log 77 = 1.97 ± 0.58 (median 1.95) towards HS 1700+6416. 
Furthermore, the distribution of 77 values is broader in compar- 
ison to the thermal fit result. Fitting a Gaussian to the log 77 dis- 
tributions shown in the middle panel of Fig.[T2]yields a FWHM 
of 0.67 + 0.12 for the temperature fit but 0.98 + 0.14 in case of 
pure turbulent broadening. An explanation for this finding may 
be the correlation of 77 with the H i column density due to the 
assumption of pure turbulent broadening which we will discuss 
in the following. If the estimate of 77 is on average too low for 
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high density absorbers, the distribution of rj values will become 
broader and shifted to lower valueS; 

One of the main conclusion of Fechner & Reimers (l2006h 
has been that the assumption of pure turbulent line broaden- 
ing may lead to the observed co rrelation between the strength 
of the H I absorption and rj (e.g. IShull et al.ll2004l) . Comparing 
the results of the fits considering the line temperature or as- 
suming turbulent line widths, we find indications that this cor- 
relation is less pronounced if thermal line broadening is taken 
into account. The average t] value for A^hi > lO'-'cm"^ ab- 
sorbers is log 77 - 2.02 + 0.55 (median 2.01) for the thermal 
fit and log 77 = 1.83 + 0.36 (median 1.87) in the turbulent 
case. If /7Hen = ^Hi is assumed, the median rj value clearly 
depends on the Hi column density of the absorber, while for 
12.0 < logA^Hi < 13.0 the median log// is 1.93, it is lower by 
0.06 dex for high column density lines (log A^H I > 13.0). On the 
other hand the mean value is independent of H i column den- 
sity if thermal broadening is considered. In this case the median 
is 2.01 in both column density ranges. 

The size of the scale of 77 variations is comparable towards 
both quasars. The scales are in the range 0.0041 < Az < 0.0263 
(mean Az - 0.0081 +0.0054, median 0.0061) coiTesponding to 
~ 3 - 4 /; ' Mpc up to 20 Mpc comoving, which is s lightly 
less than the scales found bv iFechner & ReimersI (120061) using 
the original spectrum fitting method. As also can be seen from 
the right panel of Fig. [12] we find no significant deviation for 
the distribution of scale sizes derived under the assumption of 
pure turbulent line broadening or thermal plus turbulent broad- 
ening. 

7. Summary and conclusions 

We have introduced an approach to take into account the ther- 
mal state of the IGM when analyzing the Hen Lya forest 
which has been resolved towards the quasars HE 2347-4342 
and HS 1700-1-6416. The intention of this procedure is twofold. 
At first, systematic errors due to the assumption of pure tur- 
bulent line broadening in the a nalysis of the He n Lya fo rest 
should be eliminated. Following Fech ner & Reimers (200g) the 
He ii/H I column density ratio rj may be underestimated in case 
of strong absorbers if turbulent line widths are assumed. This 
may lead to a spurious correlation of strength of the absorber 
with the 77 value. Secondly, investigating the combined Hen 
and H i Lyff forest could provide in principle an independent 
strategy to estimate the temperature and the thermal state of 
the IGM. 

The proce dure is based on the spec trum fitting method 
introduced by iFechner & Reimers (I2OO6V Instead of the ob- 
served optical data a spectrum computed from the Doppler 
profile-fitted Hi line list is fitted to the Hen data. Using 
Doppler profiles enables us to modify the temperature of the 
absorber as an additional parameter. Consequently, three pa- 
rameters are estimated, y and log Tq of the temperature-density 
relation as well as the column density ratio 77. 

Artificial data are used to test the procedure. We find 
that a minimum quality of the He 11 data, S /N > 20, is 
needed to derive reasonable constraints of IGM's thermal state. 
Furthermore, the Hi hne list may severely influence the re- 



sults. Ambiguities in the decomposition of blends can intro- 
duce systematic errors to the resulting fit and even one missing 
component may prevent the procedure from recovering the true 
temperature-density relation. 

Applying the procedure to the observed spectra of the 
quasars HE 2347-4342 and HS 1700+6416 leaves the parame- 
ters of the temperature-density relation unconstrained due to 
the low signal-to-noise ratio of the Hen data (S /N ~ 5). 
However, some of the fitted redshift ranges show indications 
of favoring either pure thermal or pure turbulent line broaden- 
ing. 

Taking into account in addition thermal broadening leads 
to different line widths for Hi and Hei with ^Heii < ^hi- 
Therefore, the He n column density resulting from the fit con- 
sidering a combination of thermal and turbulent broadening is 
systematically increased in comparison to a fit assuming the 
pure turbulent case, i.e. /^Heii - ^hi- We find a systematic 
offset of roughly 0.05 dex which is within the 1 cr uncertain- 
ties of the estimates. The lower the of the fit including the 
temperature-density relation in comparison to the;^f^ of the pure 
turbulent fit, the higher is the 77 value. Considering the com- 
bined sample of selected fit intervals towards HS 1700+6416 
and HE 2347-4342 the median log 77 = 1 .96 is obtained from 
the therm al state fit, slightly h igher than the median value 1 .93 
found by IFechner & Reimersi (2006) considering only lines 
with A^Hi < 10'^ ° cm"^, whereas the turbulent fit leads to 1.88. 

Alternatively, we estimate the thermal line width of the in- 
dividual absorption features more directly by optimizing the 
Doppler parameter of the Hen lines. In practice the values of 
the line temperature, the 77 value, and the scale of the varia- 
tion of 77 are fitted simultaneously. We find that 45 % of the 
Unes favor turbulent broadening. The inferred temperatures of 
the remaining 55 % of the components suffer from large uncer- 
tainties. Thus, it is impossible to derive a temperature-density 
relation. However, comparing the fit which includes a temper- 
ature optimization, to a fit assuming pure turbulent line broad- 
ening confirms the suspicion of Fechner & Reimers (2006) that 
neglecting thermal line widths leads to a correlation of the in- 
ferred 77 value with H i column density. While the pure turbulent 
fit results in lower 77 values for high column density absorbers 
(median log 77 = 1.87 for lines with log A^hi > 13.0 in compar- 
ison to 1.93 for fines with 12.0 < logA^Hi < 13.0), the thermal 
fit yields no correlation at all (median log 77 = 2.01 for both col- 
umn density ranges). This effect is also noticeable when com- 
paring the distribution of 77 values since a narrower distribution 
by 46 % is found if a temperature is taken into account. 

The inferred scales of 77 variations are roughly 4/1"' Mpc 
comoving. But also large scales up to Az ~ 0.03, correspond- 
ing to ~ 20 /i ' Mpc comoving, have been inferred. The derived 
scales appear to be independent of the assumption of the dom- 
inating broadening mechanism. 

Thus, we confirm the existence of a thermal component 
which results in narrower He n features in comparison to H i, 
affects the inferred 77 value. An increased Heii/Hi ratio is 
found if a thermal component is taken into account. At red- 
shifts 2.58 < z < 2.74 we find 77 ^ 100 in agreement the 
value found by Fechner et al. ('2006a') at the same redshifts to- 
wards HS 1700+6414 alone (a summary of the derived 77 val- 
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Table 1. Summary of 77 values derived with the method of direct 
temperature estimates in the redshift range 2.58 < z < 2.74. 



turbulent thermal + turbulent 

{log 77) median (log;/) median 

HE 2347 -4342 2.00 ± 0.63 1.85 2.21 ±0.82 1.88 
HS 1700 + 6416 1.97 + 0.58 1.95 2.11 ±0.32 2.03 
both sight lines 1.98 + 0.59 1.90 2.16 ±0.62 2.01 



ues using the method presented in Section |6] is given in Table 
[TJ. However, constraining the parameters of the temperature- 
density relation and thereby providing an independent estimate 
of the thermal state of the IGM is impossible with the present 
day quality of Hen Lya forest observations. In order to ob- 
tain significant constraints a signal-to-noise ratio of S /N > 20 
would be required which is 4 times higher than the S /N of the 
data available nowadays. 
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